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ABSTRACT 

We develop a maximum likelihood based method of reconstructing band powers of 
the density and velocity power spectra at each wavenumber bins from the measured 
clustering features of galaxies in redshift space, including marginalization over uncer- 
tainties inherent in the small-scale, nonlinear redshift distortion, the Fingers-of-God 
(FoG) effect. The reconstruction can be done assuming that the density and velocity 
power spectra depend on the redshift-space power spectrum having different angular 
modulations of // with fi 2n (n — 0, 1, 2) and that the model FoG effect is given as a 
multiplicative function in the redshift-space spectrum. 

By using N-body simulations and the halo catalogs, we test our method by com- 
paring the reconstructed power spectra with the spectra directly measured from the 
simulations. For the spectrum of ft or equivalently the density power spectrum P$$(k), 
our method recovers the amplitudes to a few percent accuracies up to k ~ 0.3 ft,Mpc~ 
for both dark matter and halos. For the power spectrum of ft 2 , which is equivalent 
to the density- velocity power spectrum Pse(k) in the linear regime, our method can 
recover, within the statistical errors, the input power spectrum for dark matter up 
to k ~ 0.2 /iMpc -1 and at both redshifts z = and 1, if the adequate FoG model 
being marginalized over is employed. However, for the halo spectrum that is least 
affected by the FoG effect, the reconstructed spectrum shows greater amplitudes 
than the spectrum Pse(k) inferred from the simulations over a range of wavenumbcrs 
0.05 ^ k ^ 0.3 /iMpc -1 . We argue that the disagreement is ascribed to nonlinear- 
ity effect that arises from the cross-bispectra of density and velocity perturbations. 
Using the perturbation theory and assuming Einstein gravity as in simulations, we 
derive the nonlinear correction term to the redshift-space spectrum, and find that the 
leading-order correction term is proportional to /i 2 and increases the ^t 2 -power spec- 
trum amplitudes more significantly at larger fc, at lower redshifts and for more massive 
halos. We find that adding the nonlinearity correction term to the simulation Pse(k) 
can fairly well reproduce the reconstructed Pse(k) for halos up to k ~ 0.2 hMpc^ 1 . 

Key words: cosmology: theory - galaxy clustering - dark energy - gravity test 



1 INTRODUCTION 

Cosmic accelerating expansion is the most tantalizing prob- 
lem in modern cosmology and physics. Within the frame- 
work of Einstein's general relativity (GR), the cosmic ac- 
celeration requires that roughly 70% of total energy of the 
present-day universe is in the form of unknown, mysterious 
energy component having negative pressure, dubbed as dark 
energy. An alternative explanation is the so-called modified 
gravity scenario, where the cosmic acceleration is conjec- 
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tured as a result of breakdown of Einstein's gravity on cos- 
mological scales. There are growing attempts in the com- 
munity trying to develop a consistent model of modified 
gravity that can explain the cosmic acceleration on cosmo- 
logical scales, yet recovering GR on small scales such as so- 
lar system scales, wit hout the need of dark energy (e.g. see 
iJain fc Khounfald , for a review). 

There are various methods capable of addressing the 
nature of the cosmic acceleration: type-la supernovae, clus- 
ter experiments, galaxy clustering, and weak gravitational 
lensing. These methods are sensitive to cosmic expansion 
and structure formation histories, in a complementary way, 
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over different length scales and/or different ranges of red- 
shifts. In particular, an essential approach to discriminate 
the dark energy and modified gravity scenarios is exploring 
both the cosmic expansion history and the growth rate of 
structure format ion by combining more than two different 

[Peacock et al. 1 120061 ; 



methods above (lAlbrecht et al 



2006 



2010). 



IJain fc Zhanell2008l ; iGuzik et al.l l: 

In this paper we focus on cosmological observables 
derivable from a wide-field galaxy redshift survey. A robust 
method feasible with a galaxy redshift survey is the baryon 
acoustic oscillation (BAO) experiment, which allows us to 
infer the angular diameter distance as well as the Hubble 
expans ion rate from the measured pattern of galaxy clus- 



tering (jEisenstein et alj|2005l ; ICole et aill2005l ; iBlake et all 



l2011al ). There are many ongoing and planned galaxy red- 
shift surveys aimed at achieving the BAO experiments at 
higher precisions: the Baryon Oscillatio n Spectroscopic Sur - 
vey (BOSSfl the BigBOSS project |Schlegel et all l2009h . 
the HETDEX survejQ and the Subaru Prime Focus Spec- 
trograph (PFS) project 

Adding the redshift distortion measurement can further 



improve the cosmolog 


ical power of a galaxy redshift survey 


dPeacock et al]|200ll; 


Guzzo et ai]|2008l;JWhite et al.ll2010l; 


lYamamoto et al.ll201C 


;lBlake et aLlhoilbtlSong et alKOll). 



In real space galaxy clustering is statistically isotropic in a 
statistically homogeneous and isotropic universe. However, 
in redshift space the line-of-sight component of galaxies' pe- 
culiar velocities induces an angular anisotropic modulation 
in the clustering pattern. In a structure formation scenario 
the peculiar velocities of galaxies are caused by gravitational 
attracting force in large-scale structure and the gravitational 
field can be inferred from the observed galaxy distribution 
or directly probed by weak gravitational lensing. 

More precisely there are two kinds of redshift distor- 
tion effects. One is caused by large-scale coherent velocities 
or bulk motions of halos which are associated with large- 
scale structure of large length scales, ^ 1 Mpc. This large- 
scale r edshift distor tion in the linear regime is called Kaiser 
effect (|Kaiserl [19 870. This effect amplifies clustering ampli- 
tudes of galaxies in redshift space. It is now becoming rec- 
ognized that, even at length scales of 100 /t -1 Mpc relevant 
for BAO experiments, the Kaiser effect ceases to be accu- 
rate, and nonlinearity effects needs to be included for a level 
of precision ongoing/upcoming survey can achieve. Encour- 
agingly, however, the refined, accurate modeling has been 
developed based on perturbation theor y and /or simulations 
ilScoccimarrol 12004 iMatsubaral l2008al ; iTaruva et all 120091 . 
12011 and see references therein). In this sense the Kaiser 
effect contains a cleaner cosmological information. Hence, if 
Einstein GR is a priori assumed, adding the large-scale ve- 
locity information to the BAO constraints or more generally 
the density clustering information allo ws us to significantly 
improve geometrical constraints (e.g. Alcock fc Paczvnskil 
1 19791 ; iMatsubara fc Suto|[l99rj ; iBallinger et~ai1ll996f) as well 
as cosmological parameter estimation ( e.g. [E iscnstci n et al.l 
ll999l ; lTakada" et al. 2006; lTakada|[2fJ06l ; ISaito et al.ll2008h . 

Probably more interestingly, if the density and ve- 
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locity power spectra can be reconstructed from the mea- 
sured redshift-space clustering of galaxies without assum- 
ing any gravity theory, we can now open up a window 
of exploring properties of gravity on cosmological scales 
in a model-independent way, by comparing the recon- 
structed density and velocity power spectra, because the 



ory 


dLindei 120051: IZhang et al.l 2007 


. 120081: 


Guzzo et al. 


2008 


; IWand 20081: Yamamoto et al. 


2008; 


White et al. 


200S 


iPercival & White! 200S ; Simpson & Peacock! 


201C ; 


Sons 


l2010l: Song & Kavo 201(T Yamamoto et alj 


2010; 


Reves et all l2010l: IShaoiro et all l20ld 


). For example, Ein- 
model gives us spe- 


stein gravity or a concordance ACDIv 



cific predictions on how these two spectra are related to 
each other: the two spectra have a constant overall offset 
in the amplitudes in the linear regime. Hence, if any scale- 
dependent differences in the amplitudes are found from data, 
it is a signature of failure of Einstein gravity. 

However, a viable reconstruction method needs to be 
not much influenced by uncertainties arising from small- 
scale, nonlinear redshift distortion effect due to internal 
virial motions of galaxies with in halos, the so - called Fingers- 
of-God (FoG) effect (e.g. see Ijacksonl 1 19721 : iPeacockl Il999l : 
lHamilton|[l998l :[s coccmiarr ol l2004l ). This effect causes a sig- 
nificant suppression in redshift-space clustering amplitudes 
along the line-of-sight direction. How does the small-scale 
velocity field affect the BAO-scale clustering of galaxies? 
Here is a rough estimate on the physics. Recall that virial 
velocity dispersion for massive halos of 10 15 M Q can have 
velocities of a few 10 3 km s~ . This causes a redshift mod- 
ulation given as Az ~ vu ~ 10~ 2 (in units of speed-of-light 
c — 1), which in turn causes an apparent displacement in 
the position space as Ary = Az/H(z) ~ 30 /i _1 Mpc. This 
corresponds to Fourier modes of k = 2ir/\ ~ 0.2 frMpc -1 , 
which are indeed relevant for the BAO scales. Thus the 
real-space galaxy clustering within halos at scales smaller 
than a few Mpc blows up to large scales up to ~50Mpc in 
redshift space. Since the FoG effect arises from highly non- 
linear regime and is affected by baryonic and astrophysical 
effects, it is still very challenging to have a sufficiently ac- 
curate model needed for precision cosmology. In fact, the 
FoG effect is one of the major systematic errors in galaxy 
clustering observables. 

Hence the purpose of this paper is to develop a method 
that allows us to unbiasedly reconstruct the real-space den- 
sity and velocity power spectra of large length scales from 
the measured redshift-space clustering of ga laxies, remov- 
ing t he contamination of FoG effect (also see ISong fc Kavol 
|2010| . for a similar study). This can be done by developing a 
maximum likelihood based method of reconstructing band 
powers of the real-space power spectra at each wavenumber 
bins, including marginalization over uncertainties in param- 
eters to model the FoG effect. In this method the real-space 
power spectra at each wavenumber bins are estimated such 
that the likelihood of redshift-space power spectrum is max- 
imized, assuming that the original density perturbation field 
is a Gaussian field. The real-space velocity power spectra on 
large scales include only the information on the large-scale 
redshift distortion effect, because the spectra arise from the 
density and velocity fields at physical scales corresponding 
to the wavenumbers without the FoG effect contamination. 
Our method is analogous to the cosmic microwave back- 
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groun d (CMB) power spectrum reconstruction l|Verde et al.l 
I2003T I. By using N-body simulations of 70 realizations and 
the halo catalogs, we will carefully test the method by study- 
ing whether or not the reconstructed real-space power spec- 
tra can recover the input spectra in the simulations. 

The structure of this paper is as follows. In §[2]we review 
how the redshift distortion effect due to peculiar velocities 
causes an angular modulation in redshift-space power spec- 
trum, after briefly describing how the peculiar velocity field 
is related to metric scalar perturbations. In §|3]we develop a 
maximum likelihood method of reconstructing the real-space 
density and velocity power spectra from the redshift-space 
power spectrum. After describing the N-body simulations 
and halo catalogs in § [4] we will show in § [5] the main re- 
sults of this paper; by applying the method to the N-body 
simulations and halo catalogs, we assess accuracies of recon- 
structing the real-space power spectra with the method. § [6] 
is devoted to summary and discussion. 



2 PRELIMINARIES 
2.1 Metric perturbations 

In the Newtonian gauge the perturbed Friedmann- 
Robertson- Walker metric that has scalar perturbations can 
be fully specified by the form of 



ds 



-(1 + 2*)dr + o 2 (l - 2$)dx 2 



(1) 



where a(t) is the expansion scale factor. Note that we here 
assumed a flat universe for simplicity. The metric form {TJ is 
fully general for any metric theory of gravity, as long as the 
vector and tensor perturbations are negligible. Redshift z is 
the most important observable in astronomy, and it is given 
as 1+z — l/a(t e ), where a(i e ) is the scale factor at the epoch 
when an object of interest, e.g. galaxy, emitted the photon to 
be observed by an observer. We use the convention a(to) = 1 
at present. ^ corresponds to the Newtonian potential that 
describes the acceleration of particles, while $ denotes the 
curvature perturbation. 

The expansion history of the universe is specified by the 
function of a(t) or the Hubble function H(t) = a/a, where ' 
denotes the derivative with respect to time t. Given gravity 
theory, the time evolution of a(t) or H(t) is specified once 
the energy content of the universe is specified, as in the case 
of Einstein gravity. 



2.2 The case of Einstein gravity 

Although the rest of this paper does not assume any the- 
ory of gravity, it would be instructive to discuss the case 
of Einstein gravity. This subsection also gives a background 
motivation of our work. 

Theory of gravity relates the metric perturbations in 
Eq. |T]) to matter variables. In the matter dominated era, if 
assuming the Einstein gravity, the Einstein equations yield, 
for example, the Poisson equation on sub-horizon scales, 
which relates the metric perturbation $ to the density per- 
turbation field of total matter as 

- fc 2 $ = 4ivGa 2 6. (2) 

Note that the Poisson equation here is given in the Fourier 



space, yielding the factor k 2 on the l.h.s. The matter distri- 
bution can be inferred from galaxy surveys or weak lensing 
surveys. 

In a case that the anisotropic energy stress is negligible 
as in a cold dark matter (CDM) dominated structure for- 
mation model, the two metric perturbations are equivalent 
to each other on sub-horizon scales: 



(3) 



Thus the two metric perturbations have only one degree 
of freedom, which corresponds to the density field S in the 
matter sector. 

The geodesic equation for a test particle is given by 



dX 



+ ^^P P 



0, 



(4) 



where F is the Christophel symbols. Let's consider a test 
particle which only slowly move with respect to the comov- 
ing coordinates i.e. a non-relativistic particle. Dark matter 
and galaxies are such particles. The equation of motion for 
such a test particle is given in the linear regime as 



dv I _ _ 

— + Hv = — V*, 

at a 



(5) 



where v is the comoving peculiar velocity defined as v = 
dx/dt. Thus the velocity field follows the gravitational po- 
tential. For this reason the peculiar velocity field of galaxies 
is expected as a powerful tool for probing the gravitational 
potential field. 

Another important observable is gravitational lensing. 
Solving the geodesic equation for a photon, which is a rela- 
tivistic particle, leads the lensing deflection angle to be given 
as 



Jd X Wgl(x)Vx (* + $). 



(6) 



where VFgl(x) i s f ne lensing geometrical kernel that depends 
on the backgroun d metric quantity, i.e the scale factor (e.g. 
iGuzik et al1l20l0ft . Thus lensing depends on a combination 
of the two metric perturbations, <I> -f $. 

Therefore combining different observables such as the 
galaxy distribution, the peculiar velocity and weak lensing 
in principle allow to test the consistency relation ^ = $ or 
more generally explore properties of gravity on cosmologi- 
cal scales (e.g. Ijain fc Zhanel l2008L However in this paper 
we address we can use the measured clustering features of 
galaxies in redshift space to reconstruct the power spectrum 
of the peculiar velocity field v, independently from the den- 
sity power spectrum. Hence, our method allows us to use the 
redshift-space clustering to test gravity theory on cosmolog- 
ical scales, by comparing the reconstructed power spectra of 
density and velocity fields. 



2.3 Redshift-space power spectrum 

What we can measure from a spectroscopic survey of galax- 
ies is angular positions and redshifts of the galaxies. How- 
ever, the observed redshift of a given galaxy, z, is modulated 
from the true redshift, z, due to its peculiar velocity as well 
as the metric perturbations - the so-called redshift-space 
distortion. According to the metric theory of gravity, the 
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observed redshift is given (e.g., see Sasaki 1987) as 

1 + 2 ~ (l + z){l + [*+«,]•} 
~ (l + z)[l+v x \l], 



(7) 



where <3> is the gravitational potential perturbation (Eq. [I]), 
v z denotes the line-of-sight component of the comoving pe- 
culiar velocity of tracer considered, and the notation • • • |* 
denotes the difference between quantities at observer's and 
galaxy's positions. In the second line on the r.h.s. of the 
equation above, we assumed that the effect of peculiar ve- 
locity v x , which is the order of 10~ 3 (corresponding to 
300km s _1 ) for the large-scale coherent peculiar velocities, is 
much larger than the potential amplitude $ ~ O(10 ) for 
ACDM-like cosmologies. In other words, while we will later 
focus on the density perturbations of matter or galaxies in 
large-scale structure, the effect of the metric perturbation 
$ is safely negligible compared to the density perturbations 
on relevant length scales. In addition the perturbation con- 
tributions at an observer's position (O) only contribute to 
the monopole offset (e.g. a shift in the overall normalization 
of galaxy number density at a given redshift), therefore we 
can ignore these contributions in the following. 

Via the redshift-distance relation x( z )> the apparent ra- 
dial distance to a galaxy at redshift z, x(z), is modulated 
from the true position x( z ) due to the peculiar velocity as 



X. 



X (z) + (l + z) 



dz ' 



, (1 + z) 
X(z) + u x , 



(8) 



where u x is the normalized peculiar velocity field defined as 
u x = (l + z)v x /H(z). 

The mass conservation, or the number conservation of 
galaxies, tells that the density perturbation in redshift space, 
5 S , is related to the real-space density perturbation as 



1 + 8 S 



(1 + 5) 1 + 



9u x 
dx 



1 + 5 



9u x 
&X 



+ 0(&u, u 2 ) , 



(9) 



where in the second equality of the equation above we have 
used the Taylor expansion of (1 + du/dx)~ , and we have 
ignored the higher-order terms of the perturbations (see be- 
low for further discussion). Exactly speaking the Jacobian 
transformation above breaks down when the particle mo- 
tions have shell crossing or multi-streamings at a single spa- 
tial position, which can occur in the nonlinear stage such as 
a region within a virialized halo. In other words, the equa- 
tion above is valid only at large length scales greater than 
a size of halos, which is validated on scales k 0.3 /i _1 Mpc 
we are interested in. 

Fourier-transforming the equation above yields 



5 a (k) ~ 5(k) + ^ 2 <5>(k) + o(se, e 2 ) 



(10) 



where fi is the cosine between the wavevector k and the line- 
of-sight direction. Here we have assumed that the peculiar 
velocity is irrotational, therefore is given in terms of the 
scalar velocity potential, and the quantity, 8, denotes the 
Fourier-transformed coefficient of the divergence of peculiar 



velocity field, 9 = —V • u. Also notice that in the equation 
above we have employed a distant observer approximation 
and ignored the curvature of the sky, where one axis of the 
coordinate system can be chosen to be along the line-of-sight 
direction. We again ignored the higher-order terms of the 
perturbations such as 0(58, 9 2 ). Thus the redshift distortion 
induces angle-dependent modulations, given by fi 2n (n = 
1, 2, • • • ), in the redshift-space density field. 

M otivated by the discussion abo ve and the previous 
works jKaiserll 19871 ; ISc"occimarrdl2004f ). we assume that the 
redshift-space power spectrum of galaxies (dark matter or 
halos) is given by the following functional form: 

(5 s (k)£(k')) = (2 7 r) 3 p; 5 (fc, A( )5| ) (k-k') 

PSs(k,fi) = [Pss(k) + 2f J , 2 P S e(k) + fi 4 P e e(k)] F{k,fi), 

(11) 

where Pss and Pee are the power spectra of density pertur- 
bation and velocity divergence and P$e is the cross power 
spectrum: 



<5(k)5*(k')> 
<5(k)0~*(k')) 
<fl(k)fl*(k')> 



(27r) 3 P M (fc)5| ) (k-k'), 
(27r) 3 P«,(fc),5?,(k-k'), 
(27r) 3 P ee (fc)4(k-k'). 



(12) 



The form of Eq . (Illll is often assum ed in the literature (e.g. 
lHamiltonlll998l : iTaruva et al.ll2009l . and references therein). 
However, as can be found from Eqs. (jTUJ) , pip and <|12p . we 
ignored the contributions of higher-order perturbations to 
the redshift-space power spectrum. In fact we will discuss 
later that the higher-order terms of Eq. (|10[1 can be impor- 
tant for the redshift-space power spectrum, especially for 
massive halos. To be more precise, if we recall that the ve- 
locity perturbation is smaller than the density perturbation 
at relevant low redshifts, t he leading-order corr ection to the 
Kaiser formula is found in lTaruva et al.l |2010T l to be 



5P s (k,n) 



fell U(k') 



_d 3 q_£|| 

(2tt) 3 q 2 



0(q)<5(k-q) 



(13) 



In Appendix |B] we derive the correction terms for dark mat- 
ter and halos based on the perturbation theory, and will use 
the results for the following discussion. Meanwhile we will 
assume Eq. (|11[) for simplicity. 

The function in the square bracket on the r.h.s. of 
Eq. (|lip denotes the Kaiser formula for redshift-space power 
spectrum, which is valid only at large length scales in the 
linear regime. The function F(k, fi) was introduced so as 
to take into account the nonlinear distortion effect, the so- 
called Fingers-of-God (FoG) effect, which causes a smearing 
of redshift-space clustering due to random virial motions 
of dark matter particles or galaxies within halos. Thus the 
assumption we employed in Eq. (|ll|l is the FoG redshift 
distortion and the Kaiser formula are separable functions 
in the redshift-space power spectrum. This does not neces- 
sarily hold, although an empirical model based on the halo 
model give s such a fun c tional form of redshift-space power 
spectrum (j White! l200l] ; ISeliakllioOl] . also see Hikage et al. 
in preparation). Hence the validity of Eq. (Ill[) needs to be 
further tested in combination wit h simulations. The recent 
study done in ITaruva et all (2009) gives a possible verifica- 
tion on this treatment, where it was shown that the form 
|TT} can well reproduce the simulation results in the weakly 
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nonlinear regime down to k ~ 0.2 foMpc -1 if an appropriate 
function F(k, fi) is employed. 

A theoretical understanding of the FoG effect is still 
lacking due to complicated physics involved in the nonlinear 
clustering regime. In this paper we rather employ an empiri- 
cal approach: we will consider the following functional forms 
of F(k, /i) in order to study how the results change with the 
different FoG models: 



cxp[— o 2 k 2 p?\ 
F ( k ,v) = \ l + (J 2 fcV' 



(14) 



, 2,22.^4,44 

1 — a k n + —t k fi . 



All the models have a limit of F — > 1 when k — > 0. The 
first and second forms correspond to the Gaussian and 
Lorentzian FoG mo dels that are som etimes employed in the 
literature (e.g. see iHamiltonl [l998l . for a review). We will 
treat a appearing in the forms as a free par ameter in the 
follow ing analyses, motivated by the results in iTaruv a et al.l 
(2009). The third form can be considered as a more gen- 
eral form, in analogy to the Taylor expansion of the FoG 
function in terms of k/i, and this includes the Gaussian and 
Lorentzian models in the range k/j, <C 1. Similarly we will 
treat a and r as free parameters in the model fitting. We 
will refer to these models as Gaussian, Lorentzian, Taylor-cr 
and Taylor-(a + r) models, respectively. 

Besides the assumed form of redshift-space power spec- 
trum and the FoG function (see Eqs. [11] and [14]), in the 
following we will explore a model-independent reconstruc- 
tion of the density and velocity power spectra Pss, Pse and 
Pee at each k bins, from the measured galaxy distribution 
in redshift space. More exactly speaking, since the recon- 
structed power spectra are not necessarily same as the den- 
sity and velocity spectra, our method recovers the real-space 
power spectra that are proportional to fi 2n (n = 0, 1, 2) in 
the redshift-space power spectrum, being marginalized over 
uncertainties of the FoG effect. Then we will assess the per- 
formance of this reconstruction method by comparing the 
reconstructed spectra with the spectra directly measured 
from simulations. This reconstruction problem is not a lin- 
ear problem, because the FoG function is non-linearly cou- 
pled with the density and velocity spectra. Hence, the recon- 
structed band powers at different fc-bins become correlated 
with each other even if the underlying density and velocity 
fields are Gaussian. 

Finally we remark on the Einstein gravity case, where 
the two metric perturbations are equivalent: ^ = $ as dis- 
cussed in § 12.11 In this case the density and velocity power 
spectra are related to each other in the linear regime as 
P se ~ PP SS and P ee ~ f3 2 P ss , where P = (l/b)dlnD/d]na 
with D and b being the linear growth rate and the linear 
bias parameter, respectively. Note that the possible non- 
linear correction terms (see Eq. [T3]) can be also accu- 
rately computed based on perturbation theor y and /or sim- 
ulations for a given cosmo logical model (e.g. iTaruva et al.l 
l2010l : Ijennings et alj|20lih . Thus, if the Einstein gravity is 
a priori assumed, measuring the density and velocity power 
spectra helps to break parameter degeneracies, especially the 
degeneracy between the galaxy bias and the power spectrum 
amplitudes, which in tur n helps to significan tly improve pa- 
rameter constraints fe.g. iTakada et al1l2006l ). 



3 A MAXIMUM LIKELIHOOD 

RECONSTRUCTION METHOD OF 
REDSHIFT-SPACE POWER SPECTRA 

In this section we develop a method for reconstructing the 
real-space power spectra from the galaxy distribution in red- 
shift space, based on a maximum likelihood method. 

We start with assuming that the mass density fluctua- 
tion field <5 m (x) in redshift space obeys the Gaussian likeli- 
hood function: 



£[<5 s (x)] oc 



Vdet(C) J V s J V s 



x exp 



-5 s (Xi)(C )ij<5s( X i) 



(15) 



where V s is the survey volume, C(xi — Xj) is defined as 
C(xi — Xj) = (5 s (xi)(5 s (x_ ) )), the two-point correlation func- 
tion between the density fields 5„(xi) and 5 s (xj) in redshift 
space, and C _1 is its inverse matrix. 

Analogously to t he likelihood of the CMB temperature 
power spectrum fe.g.. lVerde et al.ll2003r ). by converting the 
likelihood function to Fourier space, we can derive the log- 
likelihood function for the redshift-space power s pectrum 
(see A ppendixlAlfor the detailed derivation; also see lPercivall 
(2005J) for the similar discussion): 



•21n£= N(h,fJ,a) 



P a (h,fi a ) ln P s (k l:f i a ) _ 1 

P»(ki,fJ. a ) P s {h :f la) 



(16) 

where P s (ki, jj, a ) is the power spectrum estimated at the bin 
{ki,/j,a): 



P s {k l , fi a ) 



N(ki,n a ) 



E Nk)| 



(17) 



kG(fei,Ma) 



The quantity N(ki,[i a ) is the number of independent 
Fourier modes confined within the bin (fc;,/i a ): iV(fci,/z a ) = 

E 



1. If a surveyed volume has a cubic geome- 
try with side length L, i.e. V B = L 3 , the fundamental 
mode to discriminate different Fourier modes has the length 
given by kf = 2n/L. Hence the number of independent 
Fourier modes, for the bin (ki,^i a ), is approximately given 
as N(ki,fl a ) ~ 2nk 2 AkAij,/(2n/L) 3 in the limit h > kf, 
where Ak and A/i are the bin widths. In the equation above, 
we ignored observational effects due to survey geometry and 
masking of the surveyed region for simplicity. For actual 
data we need to include these effects. 

In Eq. (|16|l P s {ki,^ a ) is the underlying true redshift- 
space power spectrum at the bin (fei,/i a ). We assume the 
form given by Eq. pip for P s (ki,fj, a ), which is given by 
the model power spectra Pss(k), Pse(k) and Pge(k) and 
the parameters to model the FoG effect (see Eq. |14)'). 
Hence, given the measured redshift-space power spectrum 
P"(ki,fj, a ) (Eq. [17]), we can estimate the best-fit power 
spectra Pgs{k), Pse(k), and Pee(k) at each fc^-bin, includ- 
ing marginalization over the band powers at different fc-bins 
and the FoG effect parameters, in such a way that the log- 
likelihood (|16p is maximized. This is the maximum likeli- 
hood method for reconstructing the real-space spectra. 

We will demonstrate how the method above allows a 
reconstruction of the real-space power spectra using simula- 
tions. To do this, we will use the Markov-Chain Monte Carlo 
(MCMC) sampling method fe.g. lLewis fc Bridlejl2002r i. more 
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specifically Metropolis-Hasting algorithm in our work. The 
ch ain convergenc e is di agnosed by using the criteria given 
in iDunklev et all l|2005l ). The free parameters are: the band 
powers at each k bin, Pss(ki), Pse(ki) and P$e(ki), and the 
parameters to model the FoG effect given by Eq. (|14|) . where 
ki denotes the i-th wavenumber bin and the index i runs over 
the number of bins. If we employ iVbm for the bin number 
over fe m in ^ k sj fc max , the total number of model parameters 
are 3 x iVbm plus the number of the FoG parameters, 1 or 2 
(a or a and t, respectively), depending on which FoG model 
to use. In the MCMC parameter search, we adopted the fol- 
lowing priors on model parameters: Pss > Pse > Pee > 
and the FoG function < F(k,fi) ^ 1. Note that the lat- 
ter prior, < F(k, fi) ^ 1 is automatically satisfied by the 
Gaussian and Lorentzian FoG models in Eq. (|14|) . 

Another assumption we employed in the log-likelihood 
function is the Gaussian field assumption. This Gaussian 
assumption breaks down in the weakly nonlinear regime, in- 
deed over a ran ge of wavenumbe r s relev ant for the method 
above. However, iTakahashi et all (|201ll ) showed that, using 
5000 N-body simulation realizations, the non-Gaussianity of 
the density field does not cause any large impact on param- 
eter estimation in the weakly nonlinear regime. Hence we 
do not think that the non-Gaussianity affects the following 
results. 



4 N-BODY SIMULATIONS AND HALO 
CATALOGS 

To test the performance of the power spectrum reconstruc- 
tion method described in the preceding section, we will 
implement a hypothetical experiment: we will apply the 
method to mock data from N-body simulations, and then 
compare the reconstructed spectra with the input spectra 
directly measured from simulations. In this section we de- 
scribe some details of N-body simulations and the halo cat- 
alogs we will use in the following sections. 

4.1 N-body Simulations 

The N -body simulatio ns are generated by running the GAD- 
GET (|Springe]||2005h assuming a flat universe; the matter 
density f2 m = 0.238, the baryon content Ob = 0.041, the 
Hubble constant Ho = 73.2kms _1 Mpc _1 , the spectral index 
n s — 0.958 . and the amplitude of the linear power spectrum 
as = 0.76 l|Spergel et alll2007f). The trans fer function is cal- 
culated by the CAMB (|Lewis et alj|2000h . We include 512 3 

_ o 3 

N-body particles in a box of lh Gpc volume. We started 
the simulations from the initial redshift z = 30, and set the 
initial conditions of N-body particles using the Zel'dovich 
approximation. In this paper we use the outputs of z = and 
1. Our initial redshift may not be sufficiently early to accu- 
rately compute the nonlinear clu s tering of N-body particles 
as discussed in e.g. ICrocce et all (|2006h . However, the main 
purpose of this paper is to study whether the reconstruction 
method of the power spectra Pss and Pse can reproduce the 
spectra directly measured from simulations, so the accuracy 
of N-body simulations is not our concern. We will use 70 
realizations in order to reduce the statistical scatters. 

Using these N-body simulations, we also construct halo 
catalogs by adopting the friend-of-friend (FOF) method 



with the linking length of b = 0.2 (20% of the mean sep- 
aration). The minimum number of member particles is set 
to 20, which corresponds to the mass threshold of halos 
9.8 x 1O 12 M for both the z — and 1 outputs, and the re- 
sulting number density of halos is n ~ 3.8 x 10" 4 /i 3 Mpc- 3 , 
which is comparable to the number density of SDSS lumi- 
nous re d galaxies (LRGs ) targeted for the ongoing BOSS 
survey l|White et alfeOlOT) . We will use these halo catalogs to 
compute the halo power spectrum, and then address whether 
the power spectrum reconstruction method can also work for 
the halo power spectrum. 



4.2 Power Spectrum Measurement from 
Simulations 

From the simulation data above, we measure the redshift- 
space power spectrum P s (ki,n a ) in the two-dimensional 
(ki,/j, a ) bins, as well as the real-space spectra; the density- 
density power spectrum Pss(k), the density- velocity power 
spectrum Pse(k), and the velocity- velocity power spectrum 
Pee(k) (see Eq. [J2])- In the following we describe how we 
measure these spectra from the N-body simulations and the 
halo catalogs. 



4-2.1 Dark matter spectra 

First let us discuss the spectra measured from the N-body 
simulations. For redshift-space power spectrum, the distri- 
bution of N-body particles is mapped into the redshift-space 
distribution taking into account the modulation of their po- 
sitions due to the redshift distortion, where the line-of-sight 
direction is simply taken to be in the z-axis direction in each 
simulation. Note that we here adopted the distant observer 
approximation for simplicity. The power spectrum of N-body 
particle distribution is measured using the fast Fourier trans- 
form method (FFT) . In doing this we first used the "Cloud- 
in-Cell" (CIC) interpolation method for assigning N-body 
particles to the 512 3 uniformly-distributed grids in order 
to construct the grid-based density field. Then we imple- 
mented the FFT method on the density field to obtain the 
Fourier-transformed coefficients, <5 s (k). Using Eq. (|17p . we 
estimate, in each simulation realization, the redshift-space 
power spectrum, P s (ki, jj, a ), from the Fourier coefficients of 
the density field. To reduce the statistical scatters, we will 
use the averaged power spectrum of 70 realizations, and in- 
fer the lcr statistical errors from the scatters among the 70 
realizations, which correspond to the sampling variance for 
a volume of 1 /i _3 Gpc 3 . 

Fig. [1] shows the redshift-space power spectrum (color 
scales) for dark matter (N-body particles), measured from 
the simulations of z = and 1 outputs. The redshift- 
space power spectrum is compared with the real-space den- 
sity power spectrum (contours). The figure clearly shows 
redshift-space distortion effects. The Kaiser effect due to 
large-scale bulk motions increases the redshift-space power 
spectrum amplitudes along the line-of-sight direction or 
equivalently kii , stretching the iso-contours towards the 
larger ku . On the other hands, the FoG effect squashes the 
iso-contours towards the smaller k» . Comparing the left- and 
right-panels clarifies that the FoG effect is stronger at lower 
redshifts. 
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Figure 1. Power spectra measured from N-body simulations at redshifts 2 = (right panel) and z = 1 (left), respectively. The color 
scales show the redshift-space power spectrum amplitudes as a function of k± and feii, where k± and fey are wavelengths perpendicular 
and parallel to the linc-of-sight direction (which is taken as the z-axis direction in simulations). Shown is the mean power spectrum 
among the spectra of 70 realizations, each of which has a volume of 1 [h~ 1 Gpc] 3 . The anisotropic modulations of the band powers are 
due to the redshift distortion effect due to the peculiar motions of N-body particles (see text for the details) . The spectra for 2 = shows 
a stronger FoG effect: a stronger squashed feature of the iso-contours along the kit direction. For comparison, the solid contours show 
the real-space spectra, which have isotropic contours. The contours are stepped by AlogP(fc) = 0.11. 



The real-space power spectra Pss{k), Pse(k) and Pgg{k) 
are estimated from simulations as follows. The power spec- 
trum Pss{k) is just similar to the redshift-space spectrum 
as described above, but skipping a step to compute redshift 
modulation due to the peculiar velocities. For Pse(k) and 
Pee(k), we first assign the velocity components of each N- 
body particles to the 512 3 uniformly-distributed grids based 
on the CIC method, and then use the FFT method to gen- 
erate the Fourier coefficients of the velocity fields, Sj(k). 
The velocity-divergence field is computed at each Fourier 
grid as #(k) oc k • v(k). If a larger number of grids than 
512 3 (i.e. the smaller-size grid) is used, some grids may not 
contain any N-body particle, which causes an ill-behaved 
spectrum Pse(k) at small k bins. On the other hand, if we 
use a smaller number of grids than 512 3 , the CIC interpo- 
lation causes a smoothing of the velocity power spectrum 
amplitude s at large k bins we are i ntere sted in, as carefully 
studied in lPueblas fc Scoccimarrol (2009). Hence we checked 
that the 512 3 grids are rather close to an optimal choice of 
the grid number in order to avoid these artificial effects over 
a range of scales we are interested in. We again use the av- 
eraged power spectra, Pss,Pse and Pee from 70 realizations 
to reduce the statistical scatters. 

Fig. [2] shows the real-space spectra of dark matter (N- 
body particles): Pss(k), Pse(k) and P e g(k), for the two red- 
shift outputs of z = and 1. One can clearly find the relation 



Pss > Pse > Pee, and therefore the approximation given by 
Eq. (|10|l is considered valid. 

4-2.2 Halo spectra 

Now let us move on to discussion on the the power spectra 
measured from the halo catalogs. First we need to define the 
spatial position and the velocity for each halo in the simu- 
lation. We use the center-of-mass position, computed from 
N-body particles contained within each halo, as the spa- 
tial position of the halo, while we assign the mean of mem- 
ber N-body particles' velocities to the velocity of the halo. 
Then, to get the density field for the discrete halo distri- 
bution in each realization, we adopt the Nearest-Grid-Point 
(NGP) method to assign the density field in 512 3 uniformly- 
distributed grids, after the redshift modulation due to the 
halo velocity field in redshift space are taken into account. 
Similarly to the cases for N-body particles, we computed the 
density power spectra in redshift- and real-space, P s (k,[i) 
and Pss (k) . 

On the other hand, however, the velocity related power 
spectra for the halo distribution, Pse and Pgg , require some 
caution, because it is not straightforward to define the 
continuously-varying velocity field from the halo distribu- 
tion that has a much smaller number density (typically 
~ 10" 4 [fe _1 Mpc]" 3 ) than that of N-body particles. We tried 
several interpolation methods such as the CIC and the De- 
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Figure 2. The density-density (Pss), density-velocity (Pse) an d 
velocity- velocity (Pee) power spectra at z = and 1, respec- 
tively, for N-body simulation particles. Similarly to the previous 
plot, shown is the mean spectra of 70 realizations (see text for 
the details). The statistical scatters around the mean spectra are 
sufficiently small, so we do not show the scatters here (in other 
words, the average spectra are well-converged). 



launay triangulation interpolation method for which we used 
the publicly available code from the Computational Geomet- 
rical Algorithms Library (CGAL 



http:/ /www. cgal.org/ ). 



However, we could not find a reliable result for the velocity 
power spectra at scales of interest in such a way that the 
power spectra obtained become insensitive to the interpola- 
tion method. Hence, instead of pursuing a more appropriate 
method to obtain the halo velocity field, we use the grid- 
based velocity field of N-body particles, assuming no veloc- 
ity bias between the halo and N-body particle (dark matter) 
distribution. Thus we computed the power spectra, Pse and 
Pee, combining the halo density field and the N-body par- 
ticle velocity field. We will again use the mean halo power 
spectra from 70 realizations. 

Fig. [3] shows the redshift-space and real-space spectra 
for the halo distribution, measured from the 70 realizations 
according to the method we described above. Compared to 
Fig. [U the redshift-space power spectrum of halos shows 
almost no FoG effect, because the halo spectrum does not 
include contributions from the virial motions of particles 
within each halo, and rather includes only the contribution 
from the bulk motion of each halo. 



0.50 ^ 




0.05 



0.10 0.20 
k[h Mpc" 1 ] 



0.50 



Figure 3. Top panel: The redshift-space power spectra of halos, 
similarly to Fig. [T] We used halo catalogs containing halos with 
masses greater than 10 13 /i — 1 Mq, roughly corresponding to ha- 
los hosting LRGs. Slightly jaggy contours in the plot are due to 
the smaller number of halos compared to the case of N-body parti- 
cles. Compared to Fig. [T] the halo spectra shows a less FoG effect 
than in N-body particles, because halos have only bulk- velocity 
contributions in large-scale structure. Lower panel: Similarly to 
Fig. [2] the real-space power spectra for the halo distribution. For 
the density-density power spectrum Pgg, the shot noise contami- 
nation P Bn = 1/njjalo is subtracted. Note that the density power 
spectra (solid curves) have greater amplitudes at z = 1 than at 
z = because of the greater halo biases, where we used the same 
mass threshold of M m i n Ri 10 13 h — 1 Mq for both the two redshift 
outputs. For the velocity-related spectra we used the velocity field 
defined from N-body particles, instead of the velocity field of ha- 
los, because we found the difficulty of defining the velocity field 
for halos that are too sparsely sampled (see text for the details). 
Hence the velocity-velocity spectra Peg shown here are same as 
those in Fig. [2] 



5 RESULTS 

5.1 Reconstruction of matter power spectra 

We first assess performance of the maximum likelihood re- 
construction method developed in §[3]for matter spectra, by 



using N-body simulations. We stress here again that the real- 
space density and velocity spectra used to compare with the 
reconstructed power spectra, shown in figures of this and 
following sections, are the spectra directly measured from 
the simulations, and therefore include nonlinearity effects 
arising from nonlinear clustering in structure formation. 
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Figure 6. Comparing the best-fit redshift-space power spectrum, based on the maximum likelihood method, with the redshift-space 
spectrum directly measured from simulations. The solid histograms show slices of the redshift-space power spectrum amplitudes as a 
function of fi, for a fixed k: k = 0.08, 0.16 and 0.29 ZiMpc -1 from the left to right panels, respectively. Note fj, denotes the cosine angle 
between the wavevector k and the line-of-sight direction: /.t = cos(k-kii). The best-fit spectra, denoted by the different symbols, are 
computed by inserting the best-fit band powers of Pss, Psg an d Pgg at each k bin and the best-fit FoG parameters into Eq. 1111(1 . which 
are shown in Fig. [4] For comparison, the dashed curves in each panel show the spectra computed by inserting the simulation-measured 
spectra Pss, Psg and Pgg in Eq. Illl|l . but ignoring the FoG effect, i.e. setting F(k,fi) = 1. Hence the differences between the dashed 
curves and, for example, the histograms are due to the FoG effect. The reconstructed power spectra obtained using the Taylor (a + r) 
FoG model (denoted by the cross symbols) are found to well reproduce the redshift-space power spectrum. 



To apply our method to N-body simulations, we need to 
compute the likelihood function, given by Eq. (|16p . for the 
redshift-space density field. More precisely, in Eq. (|16[) , we 
need to specify a survey volume V s , which determines the 
statistical uncertainties, and need to compute the redshift- 
space power spectrum P a (ki,fi a ) at each k- and /i-bins from 
the simulations. In the following we assume Vs = lh~ 3 Gpc 3 
and use the spectrum P s (ki,fi a ) averaged from 70 realiza- 
tions to reduce the statistical scatters for illustrative pur- 
pose. For the fc-binning, we mainly use 19 wavenumber bins 
over 0.034 < k < 0.3 /iMpc -1 . We determined the bin 
widths of k and /i such that the area Afc x kA[i in the two- 
dimensional Fourier space of (fc, fj,) is kept about constant. 
Therefore, instead of using the constant bin width, Ak is 
taken to be large at small k and gradually become smaller 
at larger k. Since we use the bin width Afj, = 0.067, we use 
Afc ~ 0.02 ZiMpc" 1 at small k bins, while we use Ak = 0.01 
around k ~ 0.2 Mpc^ 1 . We will show below the recon- 
struction results for the density-density and density- velocity 
power spectra, Pss and Psg , but not for the velocity- velocity 
spectrum Pgg , because the reconstruction of Pgg is very noisy 
due to the smaller amplitudes, i.e. the small signal-to-noise 
ratio s, compared to Pss and Pgg (also see iTegmark et al.l 
|2004| . for the similar discussion). 

Fig. |3] shows the results when using the simulation 
outputs at z = 0, for dark matter (N-body particle) dis- 
tribution. The top-dotted curve shows the input density- 
density power spectrum, Pss(k), directly measured from N- 
body simulations (the average of 70 realizations). The three 
symbols around the curve, although almost perfectly over- 
lapped with each other, show the reconstructed power spec- 
tra assuming different FoG models (Eq. [14]). Note that we 
here show the results for the Gaussian and Taylor-type FoG 
models, and do not and will not show the results for the 
Lorentzian FoG model for illustrative purpose. The results 
for the Lorentzian FoG model is very similar to the results of 



the Gaussian and Taylor (a) models; that is, we have found 
that all the results for one-parameter FoG models are simi- 
lar. The error bars around the symbols, although again over- 
lapped, show la statistical uncertainties in the band power 
reconstruction at each fc-bin, including marginalization over 
uncertainties in the reconstructed band powers at different 
k bins and for different spectra (Pss, Psg, Pgg) as well as 
the FoG effect parameters. Encouragingly our method can 
well recover Pss(k), rather irrespectively of the assumed FoG 
model. 

To be more precise the upper panel of Fig. [5] shows 
fractional differences between the input and reconstructed 
spectra: AP/P = [P(reconst.) — P(input)] /P(sim.), where 
P(reconst.) and P(input) are the reconstructed spectrum 
and the spectrum directly measured from simulations, re- 
spectively. Our method recovers the input power spectrum 
within the statistical errors, achieving a few percent accu- 
racy up to k ~ 0.2/iMpc -1 . If we employ the Taylor-type 
FoG model including the orders up to (kfi) 4, (hereafter Tay- 
lor "a + r model") in Eq. (|14[) . our method can recover Pss 
up to k ~ 0.3/iMpc" 1 , which is well in the nonlinear regime. 

The lower curves with different symbols in Fig. 2] show 
the reconstruction results for the density-velocity power 
spectrum, Psg(k), assuming different FoG models (Eq. |14|). 
The reconstruction of Psg i s noisier than in Pss due to the 
lower signal-to-noise ratios (jTegmark et alj|2004 ). Also the 
reconstruction is sensitive to which FoG model is assumed, 
reflecting that the redshift-space power spectrum is affected 
by the FoG effect over a range of wavenumbers we consider. 
Fig. [4] shows that the reconstructed Psg is in closest agree- 
ment with the input spectrum, if using the Taylor-(a + r) 
FoG model that is given by two free parameters and has 
more degrees of freedom to describe a scale-dependent FoG 
effect than other models (that respectively has only one free 
parameter) . 

The lower panel of Fig. [5] shows fractional differences 
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Figure 4. The symbols show the reconstructed power spectra, 
Pgs(k) and Pgg(k), based on the maximum likelihood method 
(see Eq. 1161 and §0, for dark matter. The different symbols are 
the results assuming the different models of FoG effect: F(k, fi) 
(Eq. |14j ) in the model redshift-space power spectrum (Eq. [11]). 
The triangle symbols are the results ignoring the FoG effect F = 
1; the square shows the results assuming the Gaussian F(k,fi), 
which has a single parameter; the diamond and plus symbols show 
the results assuming the Taylor expansion forms for F(k, fi) up to 
different orders of (kfi), which are characterized by one (a) and 
two (cr, r) free parameters, respectively. The reconstructed band 
powers at each k bins include marginalization over uncertainties 
in reconstructing band powers of Pss, Pse and Peg at different 
k bins as well as the FoG parameter(s). The error bars around 
the symbols denote statistical uncertainties in the reconstruction 
for the volume 1 h~ 3 Gpc 3 , computed from the MCMC-based 
posterior distributions. Note that the band powers at different 
k bins are correlated. For comparison, the dotted curves show 
the power spectra, Pss(k) and Psg(k), directly measured from 
the simulations. For Pgg(k) the dotted curve and the symbols are 
almost perfectly overlaid. 



between the input and reconstructed spectra for Psg. Com- 
bined with the result for Pss shown in the upper panel, one 
can notice that, although Pss and Pse are unbiasedly recov- 
ered regardless of the FoG models at small k, the results 
are substantially different at large k depending on which 
FoG model to use. The FoG redshift distortion increasingly 
affects the power spectrum with increasing k. As a result, 
the reconstructed band powers at different k bins are cor- 
related with each other via the FoG effect being marginal- 
ized over, and therefore the correlations need to be properly 
taken into account (see below). Fig.[5]shows that, among the 
different FoG models, the performance of the Taylor (a + r) 
model is of promise; it can unbiasedly recover Pss over all 
the scales as well as Pse within the statistical uncertainties 
up to k ~ 0.25 ftMpc -1 , implying that the FoG model can 
nicely fit the strong FoG effect in simulations. 

To have more insights on the results in Figs. [4] and [5] 
Fig. [6] show slices of the redshift-space power spectrum am- 
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Figure 5. Using the results in Fig. [4] the plot shows frac- 
tional differences between the reconstructed power spectra and 
the spectra measured from the simulations for P$s (upper panel) 
and Psg (lower), respectively. To be more precise, AP/P = 
[P(reconst.) - P(input)] /P(input), where P input and P re const. 
are the input and reconstructed power spectra, respectively. The 
different symbols are as in the previous plot. 



plitudes, P a (k,jj,), as a function of the azimuthal angle fj,, 
for a fixed radius k: k = 0.08, 0.16 and 0.29 /iMpc -1 from 
the left to right panels, respectively. The histograms in each 
panel show the band powers measured from simulations, 
which are compared with the best-fit power spectra (dif- 
ferent symbols) obtained in Fig. [4] The best-fit spectra are 
computed by inserting the best-fit parameters (band powers 
at each k bins and the FoG parameters) into Eq. JTT} . As in 
Fig. [4] the different symbols are computed for the different 
FoG models. For comparison, we also plot the spectra, by 
dashed curves, which are computed by inserting the directly- 
measured Pss, Pse and Peg in Eq. (|llfl . without FoG term, 
i.e. F(k, j-i) = 1. Hence, the difference between the dashed 
curve and the solid histogram clarifies how strongly the FoG 
affects the redshift-space power spectrum at each k bin. 

The two extreme cases are very distinctive. At small 
k (k — 0.08ft" 1 Mpc) the FoG effect is very small, and all 
the reconstructed power spectra can well match the input 
spectra independently of the FoG models. One the other 
hand, at the largest k (k = 0.29h~ Mpc), one can clearly 
see that the FoG effect is so strong that none of the Tay- 
lor (it) or Gaussian models (also or Lorentzian model) can 
fit the ^-dependence of redshift-space power spectrum. It is 
essential to add an additional parameter in the FoG model, 
like Taylor (cr + r) model, to reproduce the simulation re- 
sults. The middle panel shows the result at the intermediate 
scale (k — 0.16/i _1 Mpc), where the FoG effect is mild and 
the FoG models of one parameter work to a good approxi- 
mation. 

Fig.[7]shows the posterior, marginalized distributions of 
the band powers, Pss, Pse, and Pee at k = 0.16 /i _1 Mpc and 
the cr parameter of the Gaussian FoG model, which are ob- 
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Figure 7. The posterior distributions of the band powers Pss, 
Pse and Pgg at k = O.I6/1 Mpc and the FoG parameter cr, for the 
power spectrum reconstruction of z = using the Gaussian FoG 
model as a demonstration example. The histograms are computed 
from MCMC chains. The solid curve in each panel represents 
a Gaussian distribution with the mean and variance given by 
the MCMC chains. The distribution of Pgg includes a range of 
Pgg = 0, meaning that the band power is not well constrained. 



tained from the MCMC chains. Here we show the reconstruc- 
tion results assuming the Gaussian FoG model, which well 
works at the scale of k = 0.16/i~ 1 Mpc as shown in Fig. [6] 
The figure shows that, while the distribution of Pss looks 
Gaussian, the distributions of Pse, Pee and a show skewed, 
non-Gaussian distributions. In particular, the distribution 
of Pgg has a wide distribution and includes a region around 
Pgg = 0, showing no constraint on the band power of Pgg. 
Given these results, we conclude that it is very difficult to re- 
liably reconstruct Pgg based on the method developed in this 
paper, at least for a survey with survey volume ~ lGpc 3 . 

The origin of the skewed distributions in Fig. is ex- 
plored in Fig. [8] which shows the posterior distributions in a 
two-parameter sub-space between the parameters in Fig. 
The figure clearly shows that the different parameters are 
correlated with each other after the nonlinear reconstruc- 
tion. In particular, the a parameter of Gaussian FoG model, 
shown here as an example, shows a strong correlation with 
the band power Psg . Thus the band powers of different spec- 
tra (Pss, Pse, Pee) at different k bins are correlated with each 
other, and the correlation needs to be properly taken into 
account for the power spectrum reconstruction. 

Given such strong correlations between the band pow- 
ers, how sensitive is the reconstruction of the power spec- 
trum to a choice of the maximum wavenumber fc max ? Fig. [9] 
studies this question. With increasing fc max , the redshift- 
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Figure 8. The color scales represent the marginalized 2d proba- 
bility distribution between the band powers (Pss, Pse, Pee) ar >d 
a as in the previous figure. The two contours levels represent the 
confidence levels of 68% and 95%, respectively. 



space power spectrum to use for the reconstruction is more 
affected by the FoG effect, and in turn the reconstructed Pss 
and Pse become increasingly affected by the FoG effect af- 
ter marginalization. The figure shows the reconstructed Pss 
and Pse obtained when including the redshift-space power 
spectrum information up to fc max = 0.16 and 0.3 ZiMpc -1 , 
respectively. For the triangle and square symbols, we as- 
sumed the Gaussian FoG model for comparison. The re- 
sults for Pss agree and are only slightly different at scales 
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Figure 9. Upper panel: Sensitivity of the reconstructed power 
spectra P$s and Pgg at z = to the maximum wavenumber 
fcmax, where the redshift-space power spectrum information up 
to fcmax is used for the power spectrum reconstruction. The tri- 
angle and square symbols show the results for fc ma x = 0.16 and 
0.30 ft _1 Mpc, respectively, assuming the Gaussian FoG effect as 
a working example. From comparison, the star symbols show the 
results assuming fcmax 

= 0.16 h~ 1 Mpc and the Taylor (a) FoG 
model. Middle and lower panels: The fractional differences be- 
tween the input and reconstructed power spectra as in Fig. \E\ 



FoG model 


fcmax = 0.3fe- 1 MpC 


fcmax = O.Wh^MpC 


Gaussian (cr) 
Taylor (a) 
Taylor (a + r) 


1=^+24.5 
1U0 34 4 

(287_ 17 ' , 356_ 76 ; 2 ) 


zi0 -137 
11 "-109 



Table 1. The best-fit FoG parameters assuming different fc m ax for 
the results of z = simulations. The units of the numbers shown 
here are /iMpc - 1 . For the Taylor (cr + r) model the parameters cr 
and t are not well constrained if using the redshift-space spectrum 
information up to fc ma x = 0.16 h~ 1 WLpc, and therefor are not 
shown here. 



around fc max = 0.16 /iMpc - , while the results for Pse are 
systematically different. Since the Gaussian FoG model can- 
not well describe the FoG effect seen in simulations as im- 
plied in Fig. IH the inaccuracy of the Gaussian FoG model 
causes a systematic underestimation in the band powers of 
Pse if including the higher-fc modes. However, the two re- 
sults for Pss and Pse agree over an overlapping range of 
fc, up to k = 0.16 /iMpc -1 , within the error bars. For 
comparison, we also show the results obtained assuming 
fcmax = 0.16 /iMpc -1 and the Taylor (cr) FoG model (see 
Eq. [14]). The Taylor (cr) model is found to give a less bi- 
ased reconstruction of Pse, implying the importance of the 
assumed FoG model even around k ~ 0.16 ftMpc -1 . 

Table [T] summarizes the best-fit FoG parameters and 




Figure 10. The contours represents the marginalized 2d prob- 
ability distribution in the parameter space (cr, r) of the Taylor 
(cr + t) FoG model (for the results shown in Figs. [4]and[5]l. The 
two contours levels represent the confidence levels 68% and 95%, 
respectively. 



FoG model 


fcmax 


= 0.3fc _1 Mpc 


fcmax = 0.16/1" ^pc 


Gaussian (cr) 
Taylor (<r) 
Taylor (cr + r) 


17.4j 
18.6j 
(240j 


18.5 
12.2 
18.8 

-18'. 4 „o R +12.6n 
-20.2' •J'J D -14.4^ 


82.6tr 7 .4 
82.4+^ 



Table 2. The best fit FoG parameters assuming different h- 
for z = 1, as in Table [T] 



the marginalized confidence ranges that are obtained for the 
reconstructions at z = assuming different fc m a X : fcmax = 0.3 
and 0.16/i" 1 Mpc, respectively. For the Taylor (cr + r) FoG 
model, the error of r parameter is smaller than that of cr, 
because the r parameter has a stronger dependence on fc/i as 
(rfc/i) 4 , than the cr-term does. Fig. llOl shows the 2D posterior 
distribution in the (cr,r) sub-space. 

Nonlincarities in matter clustering are less significant 
at higher redshifts. Hence the likelihood reconstruction of 
power spectrum we are studying may work better for higher 
redshifts. Fig. [TT] shows the results using simulations at 
2=1, similarly to Figs. and [5] Note that we assumed 
fcmax = 0.3 /iMpc -1 as done in the z — reconstruction. 
In fact the FoG effect is smaller at z, = 1 than z = 0, e.g. 
as seen from the bottom-right panel of Fig. QT] compared to 
Fig. [6] However, the Taylor (a + r) FoG model seems still 
needed in order to better recover the simulation spectra. The 
bottom-right panel clearly shows that, although the FoG ef- 
fect is smaller around fi = compared to the 2 = result 
(Fig.(6]), the Taylor (cr+r) model better captures the simula- 
tion results around fj, = ±1, showing a stronger dependence 
of fi than the Gaussian or Taylor (cr) (or Lorentzian) models 
do. This implies that it is important to properly include the 
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Figure 11. Same as in Figs. [4] \5\ and \E\ but for redshift 2 = 1. 



scale-dependent FoG effect for the power spectrum recon- 
struction, at least up to z ~ 1 we have studied. As given in 
Table [2] a non-zero r parameter is favored to capture the 
FoG effect seen in simulations. Fig. [12] shows the 2D poste- 
rior distribution in the (<t,t) sub-space, displaying a strong 
correlation between the two parameters. 



5.2 Reconstruction of halo power spectra 

Now let us move on to the reconstruction of halo power 
spectra, which are more relevant for a galaxy survey, using 
the halo catalogs constructed from 70 simulation realizations 
(see §[4] for details). The halo clustering in redshift space is 
least affected by the FoG effect, because the halos are treated 
as points and the redshift distortion effect on halo clustering 
is caused only by their bulk motions in large-scale structure, 
not by the internal virial motion within one halo. Therefore 
we can naively expect a more accurate reconstruction of the 
density and velocity power spectra for halos based on the 
maximum likelihood method we have developed in this pa- 
per. However, unfortunately, this is not that simple as shown 
below. 

For halo power spectrum we need to take into account 
the effect of shot noise arising from an imperfect sampling of 
the density fluctuation field due to the finite number of ha- 



los. In this case the maximum likelihood for power spectrum 
reconstruction needs to be modified as 



21n£ = ^7v(fc llMa 



+ ln- 



P s (hi , fl a ) -Psn 



(18) 



where P sn = 1/n is the shot noise contamination, and n is 
the mean number density of halos. In the following we simply 
assume that the shot noise is not a free parameter and given 
by the mean number density of halos we use for the power 
spectrum reconstruction: P sn = 1/n fsee ISeliak et al.ll2009l . 
for a promising method to further suppress the shot noise 
contamination). The shot noise expression is not accurate 
for th e actually measured power spectrum (e.g. ISmith et al.1 
2007), but the residual shot noise, even if exists, primarily 
contaminates to the spectrum that is proportional to fx , i.e. 
the density-density spectrum Pss- 

However the main obstacle we have faced is that we can- 
not reliably measure the velocity field of halos (see § l4.2.2l for 
details) and cannot therefore have the velocity-related power 
spectra, which are needed to assess the performance of our 
reconstruction method by comparing with the reconstructed 
spectra Pss and Pse ■ Rather we decided to use the dark mat- 
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Figure 12. The posterior distribution in the (<r, r) parameter 
space for the Taylor (a + r) FoG model at z = 1, as in Fig. [Sj 



ter (N-body particles) velocity field instead of estimating the 
halo velocity field, assuming that the halo bulk- velocity field 
is unbiased from the matter velocity field, which is often as- 
sumed in the literature. 

Hence, before going to the halo spectrum reconstruc- 
tion, we make a simple test to study whether or not the 
power spectrum reconstruction is affected by the shot noise 
contamination. This test can be done by applying our re- 
construction method to the catalogs with reduced N-body 
particles. To be more precise we randomly select N-body 
particles from each simulation realization [z = 0) until the 
number density of particles selected becomes the same to 
the density of halo catalogs, n ~ 3.8 x 10~ 4 /i 3 Mpc -3 . Then 
by using the reduced N-body particles in each simulation 
we compute the redshift-space power spectrum taking into 
account redshift modulation due to the velocity field of each 
particle. These procedures preserve the underlying spectra 
of Pss and Pse- Thus we can compare the spectra with the 
spectra reconstructed by applying the maximum likelihood 
method to the redshift-space spectrum of reduced N-body 
particles, where the shot noise contamination is subtracted 
from the measured spectrum according to Eq. (|18[) . 

Fig. [13] shows the reconstruction results (symbols in 
each panels) for the catalogs of reduced N-body particles. 
The directly measured spectra in the left panel are similar 
to the curves in Fig. [4] although we found a small difference 
in the directly measured Pss(k ) at k Si 0.2 feMpc -1 due to 
the residual shot noise effect (|Smith et al.ll2007T l. Fig. [TJ 
clearly shows that, even in the presence of shot noise, our 
reconstruction method recovers the power spectra to a sim- 
ilar precision to the results in Figs. [3] [S] and [5] Hence we 
conclude that the shot noise is not a serious source of sys- 
tematics for our method. 

Now we move to the reconstruction of halo power spec- 
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Figure 15. Ratios of the reconstructed spectra of (j? to the di- 
rectly measured P$g(k) (as in Fig. 1141 1. for different halo catalogs. 
The two results are slightly shifted in a horizontal direction for 
illustrative purpose. The two different halo catalogs are defined 
from halos with masses greater than 9.8 X 10 12 k~'Mg (trian- 
gles) and 1.86 X 10 13 h~ 1 Mq (squares), which contain at least 
20 and 38 N-body particles as members, respectively. For more 
massive halos, the reconstructed spectrum of /x 2 shows greater 
amplitudes. The solid and dashed curves show the results includ- 
ing the nonlinearity correction term to Pse(k) as done in Fig. 1141 
Pse(k) + SP '2 (fe)- To compute the nonlinearity corrections, we 
assumed fei = 1.6 and 1.84 for the less and more massive halo 
catalogs, respectively, but assumed 62 = for both the catalogs. 



tra. Fig [14] shows the results for halo catalogs at z = 0. 
First of all, the reconstruction can successfully recover the 
density power spectrum Pss(k) over a range of wavenum- 
bers we consider, as a result of properly correcting for the 
shot noise and the redshift distortion. The accurate recon- 
struction of Pss is relevant for the BAO experiments, and 
the results imply that our method may allow us to further 
use the broad-band shape of Pss (k) to improve cosmological 
constraints. However, in contrast to the results for N-body 
particles shown in Figs. I4land ll3l the reconstruction fails to 
recover the density- velocity power spectrum Pse (k) . The re- 
constructed Pse for halos gives higher amplitudes than the 
directly measured power spectrum, irrespective of the dif- 
ferent FoG models. In fact, as explicitly shown in the right 
panel, the measured redshift-space power spectrum shows 
greater amplitudes than predicted by the Kaiser formula 
(Eq. [TT] with no FoG effect, i.e. F = 1). The enhancement 
in the power spectrum amplitudes is opposite to the FoG 
effect, which always suppresses the amplitudes. 

We argue below that the results in Fig.[l4]can be under- 
stood by the nonlinearity effect on the redshift-space power 
spectrum. As we briefly discussed around Eq. (JT3} , the non- 
linear clustering causes a correction to th e Kaiser formula of 
redshift-space power spectrum (also see IScoccimarrol |2004| ; 
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Figure 13. Reconstruction results of the power spectra for the reduced N-body particle distribution as in Fig. 1111 where a smaller 
number of N-body particles are randomly selected in each simulation realization in such a way that the mean number density of the 
resulting particles becomes comparable to that of halo catalogs we will use below (see §gj): n= 3.8 x 10" 4 h 3 Mpc- 3 . In this case the 
shot noise term arising from a finite number of the sampled N-body particles affects the redshift-space power spectrum measurement. We 
applied the power spectrum reconstruction method (Eq. 1161 ) to the redshift-space power spectrum after simply subtracting the expected 
shot noise term 1/n from the measured spectrum. The dotted curves in the left and right panels and the spectrum in denominator in 
the middle panel are the input power spectrum, which are the same as the spectra of original N-body particles in Figs. [4] [5] and [6] It is 
found that our method nicely recovers the input power spectra even in the presence of shot noise contamination. 



iTaruva et aHl20ld . for a more extensive discussion). Assum- 
ing that the density perturbation is greater than the velocity 
field, which can be even more validated for highly biased ha- 
los with 6 > 1, the leading-order correction term is found to 
arise from the cross-bispectrum {586) (see Eq. |13|): 



5P s (fc,/i) <- fci 



d 3 q 



(2tt) 3 q 2 



g|l 0(q)5(k- 



(19) 



We tried to measure this correction term from the simula- 
tions, but could not obtain the reliable results as the bis- 
pectrum measurement is very noisy. Instead we here use the 
perturbation theory prediction assuming a ACDM cosmol- 
ogy (or equivalently Einstein gravity). In Appendix [Bl we 
explicitly derive the leading-order correction term given as 
a function of the linear mass power spectrum. We find that 
the leading-order correction only contributes to the redshift- 
space power spectrum at the power of fi 2 : 

PLio(k,fi) = Pss(k) + 2 fJ : 2 [P sg (k) + SP^ik)] +n 4 P ge (k), 

(20) 

where 8P^i (k) is the correction term. Including the halo bias 
parameters, the correction term is expressed as 



SP^(k) 



fbjk 3 

(27T)3 



dr 



dx x 



7 v 



7rx ■ 



6x 2 )P s L s (k) 



+-(7x + 3r 



10rx 2 )P s L s (kr) 



P s L s(Wl + r 2 -2rx) 



(1 + r 2 - 2rx) 



+ 



-ifff(fc) 



(2n) 2 



dr-(l + r 2 )P s L s (kr) 



dr / dx xrP^ s {kr)Pg S {k\/ 1 + r 2 



2rx), 
(21) 



where / = d\n D/dlna (D is the growth rate) and fei and 62 
are the linear and nonlinear bias parameters. Note that we 



here assumed that the velocity field of halos is unbiased with 
respect to the velocity field of dark matter. Eq. (|21|) clearly 
shows that the nonlinearity correction term depends on halo 
bias parameters. The first term depends on the linear bias 
parameter as oc b\. Compared to the density- velocity power 
spectrum Pse(k), which depends on 61 as Pse oc 61, the non- 
linearity correction term can be more importa nt for more bi- 
ased h alos, or equivalently more massive halos. ITaruva et al.l 
(2010) derived more comprehensive equations for dark mat- 
ter including other nonlinear terms which arise from the 
cross-bispectra such as (899). The other terms are found to 



have the contributions of [i n (n = 1, 2, ■ • ■ 4), but depend 
on halo bias as oc 61. For highly biased halos with 61 > 1 as 
halos we are studying, the term given by Eq. ()21|) has most 
dominant contribution. However, the perturbation theory is 
known to be less accurate for lower redshifts such as z = 0, 
due to the stronger nonlinear clustering effects. Hence the 
results shown here still need to be more carefully studied. 

The dashed curves in the upper-left and -right panels 
of Fig. Q3] show the results where we added the nonlinearity 
correction term to P$e(k) measured from the simulation as- 
suming the Taylor (cr + r) FoG model: Pse{k) + SP^ (k). To 
obtain the theory prediction of the nonlinearity correction 
term, we assumed fei = 1.6 and &2 = 0, where 61 is esti- 
mated by comparing the density power spectra Pss for dark 
matter and halos at small k. More exactly we computed the 
correction power sp ectrum of fj? using the full expression in 
ITaruva et alj l|2010t) . also taking into account the halo bias 
dependences on the different terms that are either propor- 
tional to 61 or b\. Eq. (|21[) gives the similar shape, but about 
10% higher amplitudes than the curves in Fig. 1141 Inter- 
estingly, the nonlinearity correction term increases, rather 
than suppresses, the amplitudes of real-space power spec- 
trum that is proportional to fi 2 in the redshift-space power 
spectrum. The enhancement increases with increasing k. We 
should also emphasize that such a nice agreement including 
the nonlinearity correction can be found only if using the 
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Figure 14. Same as in Figs. l4l l5l and l6l but for halo spectra at z = 0. As described in § 14.2.21 the input density- velocity spectra Psb, 
which are used to compare with the reconstructed spectra, we used the spectra between the halo density field and the N-body particle 
velocity field, because the halo velocity field is hard to construct due to a coarse sampling of the velocity field. As can be clearly seen 
from the lower three panels (especially two lower-right panels), the measured redshift-space power spectra show greater amplitudes than 
the redshift-space spectrum inferred from the simulation, without the FoG effect F = 1 (see text for discussion). The dashed curves in 
the upper-left and -right panels show the results obtained by adding the perturbation theory prediction of the nonlinearity correction 
term (Eq. 1211 ) to the directly measured Pgg(k). In the predictions, we assumed the halo bias parameters bi = 1.6 and 62 = 0, where the 
linear bias parameter is estimated by comparing the density power spectra (P$s) of dark matter and halos at small k. 



Taylor- (a + r) FoG model. Given the fact that the redshift- 
space power spectrum of halos is least affected by the FoG ef- 
fect, our results imply that the Taylor (u+r) model has more 
degrees of freedom than other one-parameter FoG models 
and can effectively capture higher-order contributions of /i 2n 
(n = 1, 2, ■ ■ ■ ) that ari se from nonlinearity effects as studied 
in lTaruva etail l|2010l ). 

We also studied the halo power spectrum reconstruc- 
tion using the halo catalogs constructed from z = 1 sim- 
ulations. We similarly found that the reconstructed power 
spectra of fj 2 show greater amplitudes than expected from 
the measured Pse(k). The nonlinearity correction term is 
found to similarly explain the reconstructed power spec- 
trum, in slightly less agreement, if using the reconstructed 
power spectrum obtained from the Taylor (cr+r) FoG model. 
We have also found a subtle contamination of the residual 
shot noise for the z = 1 results, and therefore we here show 
the results for z = for illustrative clarity. 

To obtain more insights on the halo power spectrum 



results, in Fig. [15] we study the reconstruction results for 
Pss using different halo catalogs where halos are selected 
with different mass thresholds. To be more precise, we 
made the new catalogs by employing higher mass threshold, 
1.86x 1O 13 /i _1 M (including more than 38 N-body member 
particles), rather than the threshold 9.8 x 10 12 h~ 1 M & (20 
particles) we have so far used. Note that the number density 
for more massive halos is ft ~ 1.9 x 10 -4 /i 3 Mpc -3 , compared 
to 3.8 x 10~ 4 /i 3 Mpc" 3 for our fiducial halo catalogs. The 
estimated bias parameter is foi = 1.84 compared to 61 = 1.6. 
Fig. \TE\ shows the reconstruction power spectrum of fi 2 for 
the different halo catalogs, compared to the density- velocity 
spectrum. The figure shows that the reconstructed spectrum 
for more massive halos has higher amplitudes than for less 
massive halos. The solid and dotted curves show the predic- 
tions obtained by adding the nonlinearity correction term 
(Eq. [21]) to the directly measured Pse{k), where we used 
in the computation the linear bias parameters above and 
assumed 62 = for simplicity. The nonlinearity correction, 
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which depends on the halo bias, fairly well reproduces the 
reconstruction results. In summary such nonlinearity correc- 
tion terms need to be included when interpreting the recon- 
structed power spectra for halos, or more generally galax- 
ies. We again emphasize that our method reconstructs band 
powers of the real-space power spectrum, which are propor- 
tional to (i 2 in the redshift-space power spectrum, rather 
than the band powers of Pse alone. 

Finally we comment on the impact of nonlinearity cor- 
rection on the reconstructions results for dark matter (N- 
body particles), which we showed in the preceding section. 
For the dark matter spectrum, which has b — 1 by defini- 
tion, the nonlinearity correction (Eq. [5T]) is smaller com- 
pared to the results of halos. However, using the pertur- 
bation theory predictions, we found that the nonlinearity 
correction is not negligible. Including the nonlinearity cor- 
rection improves agreement with the input power spectrum 
over a range of wavenumbers up to k ~ 0.2 /iMpc -1 for 
the Taylor-(cr + r) FoG results shown in the middle panel 
of Fig. [5] However, the nonlinearity correction increases the 
disagreement at the larger k. In summary we conclude that 
our reconstruction method can well recover the real-space 
power spectrum, which is proportional to in the redshift- 
space power spectrum, up to k ~ 0.2 /iMpc -1 for both dark 
matter and halos, if including the nonlinearity corrections. 



6 SUMMARY AND DISCUSSION 

In this paper we have developed a maximum likelihood 
based method of reconstructing the real-space power spec- 
tra of density and velocity fields, from the two-dimensional, 
redshift-space clustering of dark matter and halos (suppos- 
edly galaxies). This method is developed in anal ogy with the 
CMB power spectrum reconstruction method l|Verde et al.l 
1200 j ). 

By assuming the form of redshift-space power spectrum 
given by Eq. (|16l) . we developed a method of reconstruct- 
ing the band powers of Pss and Pse at each k bins, being 
marginalized over uncertainties in the band powers at dif- 
ferent k bins and the parameters to model the FoG effect, 
in such a way that the likelihood of the redshift-space power 
spectrum measured becomes maximized. One assumption 
we have employed for the method is the functional form of 
redshift-space power spectrum (Eq. [TB]), where the Kaiser 
formula and the FoG effect is given by multiplicative func- 
tions. In fact this form is expected based on the halo model 
picture l| White! l200ll ; ISeliakfl200ll . also see Hikage et al. in 
preparation). The real-space power spectra, especially at 
such large length scales (k ^ 0.3 /i -1 Mpc), contains cleaner 
cosmological information in the linear or quasi-nonlinear 
regimes, and are relatively easier to develop a sufficiently 
accurate model by using a suit of simulations and/or re- 
fined perturbation theory. Furthermore, by measuring the 
velocity-related power spectra in a model-independent way, 
we can open up a new window of testing gravity on cos- 
mological scales. That is, we can address whether or not 
the velocity field inferred is consistent with the gravity field 
inferred from the density field, because the density and ve- 
locity fields are related to each other via gravity theory. 

We have carefully tested our method by comparing the 
reconstructed real-space power spectra with the spectra di- 



rectly measured from simulations of 70 realizations, for dark 
matter (N-body particles) as well as halos. For matter power 
spectra (i.e. N-body particles), we showed our method nicely 
recovers the power spectra Pss and Pse over a range of scales 
k 0.3 /i -1 Mpc and at redshifts z = and z = 1 (see 
Figs. HI 151 and ITTT) . to accuracies within the statistical errors, 
if we use the Taylor (cr + r) FoG model (see Eq. [14]), which 
has more degrees of freedom (2 parameters) than the other 
models, Gaussian, Lorentzian and single-parameter Taylor 
FoG models. Our results imply that the FoG effect seen in 
simulations has a complex scale-dependence, and is impor- 
tant to take into account the scale dependence in order to 
obtain an unbiased reconstruction of the band powers of 
Pss and Pse at scales down to k ~ 0.3/i _1 Mpc. In other 
words, the FoG effect affects the redshift-space power spec- 
trum over a wide range of wavenumbers. Hence an inaccu- 
rate modeling of the FoG effect causes a biased estimate of 
the power spectra. It is also worth noting that the recon- 
struction causes correlations between the band powers of 
different power spectra and at different fc-bins and the FoG 
model parameters (see Figs. 151 and llip . 

For the halo power spectrum, we showed that our 
method again nicely recovers the density power spectrum 
Pss over a wide range of wavelengths, up to k ~ 0.3 hMpc~ 
(see Fig. I14|) . Such an accurate reconstruction of Pss is very 
promising, because the shape and amplitude information of 
Pss are sensitive to cosmological parameters such as the tilt 
and running index o f the primordial power spectrum and 
neutrino masses (e.g. iTakada et aill2006l ; ISaito et alj|2009l . 
2010). On a measurement side, the halo power spectrum 
can be estimated from actual gal axy redshif t surv ey, e.g. 
based on the method developed in iReid et al l l|2010h where 
galaxy pairs with small spatial separations are clipped out. 
Although the halo power spectrum is supposed to be less 
contaminated by the FoG effect, it is very important to min- 
imize the residual FoG contamination in order to extract 
unbiased cosmological information from the measured halo 
power spectrum. For example, since the FoG effect causes 
a suppression in the power spectrum amplitude, a residual 
FoG contamination would cause a bias in neutrino mass con- 
straints because the main effect of massive neutrinos is also 
the suppression on power spectrum amplitude (Hikage et 
al. in preparation). Our method can give a robust way of 
measuring the density power spectrum, minimizing the FoG 
contamination. 

For the halo velocity power spectrum Pse, we found 
some difficulty. First of all, we could not reliably recon- 
struct the velocity field of halos from simulations, due to too 
sparse sampling of halos' velocities. Hence we instead used 
the velocity field of dark matter (N-body particles) assuming 
that the large-scale bulk motions of halos are unbiased from 
the velocities of dark matter, which has been often assumed 
in the literature. Note that the real-space velocity field we 
considered here contains only the large-scale information at 
k £ 0.3 /iMpc -1 , and therefore is not affected by any virial 
motions within halos. As a result, we found that the recon- 
structed power spectrum of /i 2 systematically differs from 
Pse(k) directly measured from the simulations (Fig. I14[) . 
In fact, the measured redshift-space halo power spectrum 
shows greater amplitudes than the spectrum inferred from 
the Kaiser formula of redshift-space spectrum, without the 
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FoG effect that causes a suppression in the redshift-space 
power spectrum amplitudes. 

Therefore we argued that the halo power spectrum is 
affected by the nonlinearity effect. In Appendix [B] assum- 
ing a ACDM cosmology or Einstein gravity, we derived 
the leading-order nonlinearity correction to the Kaiser for- 
mula of redshift-space power spectrum, which arises from 
the cross-bispectrum between the density and velocity per- 
turbations. We meant by the leading-order term that the 
term appears to have the largest contribution, assuming 
that the density perturbation is greater than the velocity 
field at length scales of interest. We found that the leading- 
order contribution is proportional to fi 2 and has greater 
amplitudes for more biased halos, i.e. more massive halos. 
We showed that adding the perturbation theory prediction 
to the simulation Pse(k) better matches the reconstructed 
power spectrum (see Fig. I14[) . We also found that, by us- 
ing the different halo catalogs defined with different mass 
thresholds, the halo bias dependence of the nonlinearity cor- 
rection is seen in the reconstructed power spectra of p? (see 
Fig. US}. 

Hence a more appropriate statement for our maximum 
likelihood method is that the method can recover the real- 
space power spectra, which are proportional to fi° and p 2 , 
respectively, in the measured redshift-space power spectrum, 
including marginalization over uncertainties in the FoG ef- 
fect. In other words the reconstructed spectrum of y? is not 
necessarily the same as the density-velocity power spectrum 
Pse(k), which we have used in our comparison. The nonlin- 
earity effect on the real-space power spectrum of fj 2 needs 
to be included if we want to use the reconstructed power 
spectrum to constrain cosmological parameters as well as to 
test gravity theory. On the other hand, we found that the 
power spectrum of /i 4 is very noisy to reconstruct, for the 
ranges of wavenumbers and redshifts we have considered in 
this paper. 

Recentlv lTaruva et all j|2010h studied the redshift-space 
power spectrum including nonlinearity effects, based on the 
extended perturbation theory. They found that the non- 
linear correction terms including the higher-order terms of 
0(8 2 ) have the contributions that are proportional to p? n 
(n — 1, 2, . . . , 4) in the redshift-space power spectrum. The 
comparison of the theoretical prediction with the recon- 
structed power spectrum based on our method is very in- 
teresting, and will be studied elsewhere. 

One encouraging result is our method can unbiasedly 
recover the real-space density power spectrum Pss{k) even 
in the presence of redshift distortion effect. Given this re- 
sult our method may offer even a new means of obtaining 
geometrical constraints on the Hubble expansion rate and 
the angular diameter distances beyond the usual BAO con- 
straints. We have assumed throughout this paper that the 
underlying cosmology is known. However, this is obviously 
not true for an actual observation. In reality, we have to as- 
sume a reference cosmological model to perform the cluster- 
ing analysis of galaxies, and the assumed cosmology gener- 
ally differs from the underlying true cosmology. An imperfect 
cosmological model causes additional angular anisotropies 
in the measured redshift-space power spectrum - the so- 
called cosmological distortion. In terms of Eq. (|16p an in- 
correct cosmology leads some power of the density power 
spectrum Pss (fc) of /x° to leak into the power spectrum with 



powers higher than fi 2 in the redshift-space power spec- 
trum. Contrary, if we seek the reconstructed Pss of maxi- 
mum amplitudes with varying reference cosmological mod- 
els, we may be able to obtain the cosm ological constraints 
(also see iPadmanabhan fc White| [2008. for a similar discus- 
sion). jriTis_jneH2£d_j£oks simila r to the Alcock-Paczynsk i 
test llAlcock fc PaczvnskH 1 19791 ; iMatsubara fc Sutol Il996l ; 
Ball inger et all 1 19961 ). but our method may have practical 
advantages: our method allows us to measure the real-space 
power spectrum of /i° in a model-independent way as well 
as to derive cosmological constraints being marginalized over 
uncertainties in the FoG effect. The feasibility of this method 
is our future work and will be presented elsewhere. 

Our reconstruction method is done in the two- 
dimensional Fourier space of (k,n). In practice one may 
want to exclude the Fourier modes around fj, ~ ±1, which 
are more affected by the FoG effect. In our method it is 
straightforward to include a masking of the modes around 
ft±l; that is, the real-space power spectra are reconstructed 
by using the Fourier modes in redshift space, excluding the 
modes around p±l. We have tried several masking methods 
of n, but could not find any significant differences from the 
results shown in this paper. 

In this paper we have ignored some observational effects 
for simplicity. For example, to apply our method to actual 
data, we need to include effects such as survey window func- 
tion and the cur vature of the sky. Th ese effects have been 
well studied fe.g. iTegmark et al.ll2004r ). and would be rather 
straightforward to include, although a further careful study 
needs to be done. 
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APPENDIX A: LIKELIHOOD FUNCTION OF 
REDSHIFT-SPACE POWER SPECTRUM 

In this appendix we derive the likelihood function of redshift- 
space power spectrum. To do this, we assume that the den- 
sity fluctuation field of large-scale structure tracers (matter 
or galaxies) is Gaussian and obeys the following Gaussian 
likelihood function: 



C oc 



VdetC 



exp 



d Xi f d x 



V 



V 



(Al) 

where <5(x) is the density fluctuation field (of matter or 
galaxies), dj is the correlation matrix between the fields 
<5(xi) and 8(xj), C" 1 is the inverse matrix, and detC is the 
determinant of the matrix C. The covariance matrix or the 
correlation function, Cij, can be expressed in terms of the 
power spectrum P(k) as 

Ca = <*(x i )*(x i )> 



^^P(fc)e k(x -^ 



(A2) 



where V s is the survey volume. Given a finite-volume sur- 
vey, we introduced the dis crete Fourier transform ation of 
the density field (e.g. see iTakada fc Bridle! [2007b . where 
the fundamental Fourier mode is given by the survey size; 
kf = 27r/Z/ (V s = L z ). Here we ignored survey geometry and 
boundary effects for simplicity. Therefore the inverse of the 
covariance matrix can be given as 



E 



Vs 

P(k)' 



ik- (x 



(A3) 



This can be proved because the product of the covariance 
matrix and its inverse matrix satisfies the orthogonal rela- 
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tion, which is formally computed as 



EE P «pi7T elkX,elk '' Xfc ^- k ' 

k k' V ' 



(A4) 



where 8 K (k — k') is the Kronecker-type delta function; 
8 K (k — k') = 1 if k = k' within the bin w idth, otherwise 
5 K (k~ k') = (see lTakada fc Bridle! (|2007h V 

Using the equations above, the argument in the expo- 
nential of Eq. (| All) can be reduced to the following form in 
Fourier space: 



V. 



^EEE^ k 

v s 



Vs 



P(k') 



ik' (x,- — x,,- ) r ik -x^ 



k"Ok+k" , -k 



P(k) 

\k\ 2 
p(k)- 



(A5) 



Note that, in equation above, I assumed that \Sk\ 2 and P(k) 
have same dimensions such that the combination \Sk\ 2 / P(k) 
becomes dimension-less. 

Therefore the log-likelihood function of density fluctua- 
tion field (Eq. [Al] ) is reduced to the log- likelihood function 
of the power spectrum: 



21n£ 



E 



p(k) 



lnP(fc) 



(A6) 



where we have ignored the constant additive term. The log- 
likelihood above is rewritten in terms of the power spectrum 
estimator as 



-21n£ = N{ki,» a ) 

hi ,fj. a 



P(k iy Ha) 
P(ki,H a ) 



+ In P(ki, Ha) 



(A7) 



Here P is the power spectrum estimator 

1 



P(ki, fl a ) 



N{ki,Ha 



■ E 

k£fej,f(, 



(A8) 

is over Fourier modes k 



where the summation X^keo f ) 
satisfying the condition that the wavevector k lies in the 
bin labeled by the length and the azimuthal angle between 
the line-of-sight direction and the wavevector, (ki,(J, a ), and 
N{ki,fia) is the number of independent Fourier modes: 
N(ki,fia) = E ke(fel , Ma ) « 27rfc?AfcA M /(27r/L) 3 for h » 
2n/L. The /i a -dependence of power spectrum accounts for 
redshift distortion effect. 

Adding the constant term into the equation above such 



that the log-likelihood function can be maximized if the the- 
ory power spectrum P(ki,/j, a ) is equa l to the estimated one, 
we can arrive, as in the CMB case (| Verde et al.ll2003h . at 
the expression: 



-21n£= N(k h n a ) 



P(ki,/l a ) 



+ ln 



P{ki,Ha) P(ki,Ha) 



(A9) 

Hence we can use this log-likelihood to estimate the underly- 
ing power spectrum P(ki, fi a ) at each bin, given the observed 
power spectrum P(k i: fj, a ). 



APPENDIX B: NONLINEAR CORRECTION 
TERMS OF THE KAISER FORMULA 

In this appendix , foll owing the method developed in 
IScoccimarrol (|2004) and iTaruva etU] (|201Ch . we derive the 
leading-order correction term of higher-order perturbations 
to the Kaiser formula for the redshift-space power spectrum, 
assuming a ACDM cosmological model based on the Einstein 
gravity. 

Let us begin our discussion with recalling that the red- 
shift distortion effect on a given tracer of the large-scale 
structure is recognized as a mapping between redshift- and 
real-space positional vectors: 



S = X + U-yZ, 



(Bl) 



where s and x are the positional vectors in redshift- and 
real-spaces, respectively, u x is the line-of-sight component 
of the normalized peculiar velocity (see around Eq. [5]) and 
2 is the unit vector of the line-of-side direction in the real- 
space coordinate system. The mass or number conservation 
law gives the relation between the density perturbations in 
redshift- and real-spaces: 



S s (s) = [l + tf(x) 



1 + — * x 



(B2) 



Therefore the Fourier transform of the redshift-space 
density perturbation can be expressed in terms of the real- 
space density and velocity perturbation fields as 



* B (k) 



j3 

a s 



cfx 
d 3 x 



d 3 s 5 a (s)e lk s 



8Uy 

ox 



1 + 5- 



— ikx i j~*r 2\ 

e + 0{u ) 



S(k) + ^ 2 0(k) + k/i 



d 3 q q [{ - 



(2tt) 3 q 2 



<5(k-q)0(q) + O(0 2 ), 
(B3) 



where we have expressed the peculiar velocity field as the 
velocity-divergence field as u x — i{ku/k)6 (see around 
Eq. [10]), and [i = ku/k. On the third equality on the r.h.s. 
of the equation above, we have used the Jacobian, \dx/ds\ 
to make the integration variable change, s — > x. Here, given 
the fact that the density perturbation is greater than the 
velocity perturbation, we kept the leading-order nonlinear 
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correction term which has the order of O(50) and ignored 
the higher-order terms than O(0 2 ). 

Therefore, the nonlinear correction term to the redshift- 
space power spectrum is found to be 



5P s (fc, M ) = 2(fc M ) J±± 



d q 9|| 



B Me (-k,k-q,q), (B4) 



(2tt) 3 q 2 

where the bispectrum is defined as 

<5(ki)S(k 2 )0(k 3 )> = B Mfl (ki,ka, k 3 )(27r) 3 4(k 1 + k 2 + k 3 ). 

(B5) 

Using the perturbation th eory of structure formation (e.g. 
I Jain fc Bertschinger|[l993 ) we can express the bispectrum in 
terms of the linear power spectrum as 

S Me (ki,k 2) k 3 ) = 2/ [F 2 (k2,k3)P&(fc 2 )P 4 L 4 (fc3) 

+F 2 (k.i,-k 3 )Pss(ki)Pfs(k 3 ) + G a (ki,k a )J^(*i)if,(fa)] , 

(B6) 

where / = d\nD/d\na, 9 = fS at the linear order for our 
notation and the kernels F2 and G2 are defined as 



F 2 (k!,k 2 ) 
G 2 (ki,k 2 ) 
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(B7) 



For halo clustering, we need to further take into account 
halo bias. Here we simply assume that halo bias is determin- 
istically g iven as a function of th e underlying mass density 
field fe.g. lFrv fc Gaztanaga|[l993 ): 



(5 h (x) = Mm(x) + y5„(x), 



(B8) 



where 61 and 62 are the linear and nonlinear bias parameters. 
On the other hand we assume that the velocity field of halos 
is unbiased to the velocity field of dark matter. Then we can 
similarly compute the correction term to the redshift-space 
power spectrum of halos: 



5P s , h (k, (j,; z) 

+2/&1&2 



26? (fc/*) 
d 3 q 
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(B9) 



We can further simplify the equation (|B9|l by using 
the usual fo r mula develope d in Appendix A and B of 
iTaruva et all (|20ld ) (also see iMatsubaral (^OOSbl ) 1 ). The for- 
mula tells that the following equation holds for an arbitrary 
scalar function /(q, k): 



d 3 q q« 



(2tt) 3 q 2 



/(k,q) = /xQ(k,q), 



where /i = k z /k and 
k 2 



k 2 f 00 r 1 
Q(fc) = j^pp J dr J dx rx f( k , r , x )- 



(BIO) 



(Bll) 



Here we have introduced the integration variable transfor- 
mations as q = kr and q • k = k 2 rx. 

Therefore, by comparing Eqs. (|B9[I and (|B 10[l . we can 
find that the nonlinear correction term to the Kaiser for- 
mula of the redshift-space power spectrum is proportional 
to fi 2 . That is, the leading-order nonlinear correction term 
only contaminates to the density-velocity power spectrum 
in the Kaiser formula. If we define the correction term as 
5P s Kaiser (fc,/i) = 2/i 2 5P M 2(fc), we find that the real-space 
power spectrum SP^i (k) is given as 



«>(*) = 



fbjk 3 
(2tt) 2 



dr 



dx x 



{ r —{-l + 7rx-Qx 2 )Ph{k) 



1 1 f*j . q ln 2,„l„, P S s(kVl + r 2 -2rx) 
+ -(7x + 3r-Wrx )P w (fcr)j {1 + r ,_ 2rx) 



-Ptl{h) / dr-(l + r 2 )P s L s (kr) 
Jo -5 



+ 



fb!b 2 t 
(2n) 2 



r 00 pi | 

/ dr dx xrPss(kr)P$§(k*\/ 1 + r 2 — 2rx). 

(B12) 



Eq. (|BI2[l has several interesting implications. First of 
all, the nonlinear correction term scale with halo bias as 
8P^2 oc b\. Since the density- velocity power spectrum for 
halos scales as Pse oc 61, the correction term can be more 
important for more biased halos. Secondly, once the gravity 
theory is assumed (here the Einstein gravity), the correction 
term can be comput ed as a function of cosmological models. 
ITaruva et all £2010) further derived other correction terms 
arising from the higher-order perturbations than 0(9 2 ), and 
then showed that the nonlinear correction term can show 
remarkable nice agreement with the simulation results. In 
the main text (see Sec. 15. 2|) we use Eq. (|B12|) to explain 
the power spectrum reconstruction results for halos, where 
we found that the reconstructed power spectrum, which is 
proportional to fi 2 , shows sizable difference from the input 
Pse(k). 
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